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•FOREWORD 


This work was performed under NASA Goddard Space Flight 
Center Grant NGS-5288. The basic purpose of the research 
performed was to model subsurface media by employing stochas- 
tic techniques. A preliminary investigation indicated that 
modeling subsurface media by discrete particles having random 
position and orientation would be a viable procedure. The 
method has wide applicability since most subsurface media can 
be viewed as an aggregate of dielectric particles. In addi- 
tion, several multiple scattering methods exist which enable 
one to calcualte the electromagnetic radiation from such a 
collection of particles for various parameter regimes. 

To demonstrate the usefulness of the method, the electro- 
magnetic backscattering from a half space of discrete lossy 
dielectric scatterers was analyzed. The method of Foldy was 
employed to find an equation for the mean field in the scat- 
tering region. From this equation, an effective permittivity 
was obtained. The effective permittivity was anisotropic re- 
flecting the non-spherical nature of the particles being con- 
sidered. Following this, the correlation of the scattered 
field was found by employing the distorted Born approximation. 
This method treats the scatterers as particles embedded in the 
effective medium. The backscattering coefficients were then 
computed . 

Specialized results were obtained by letting the arbi- 
trary particles take the form of discs. Since lossy dielec- 
tric discs closely resemble leaves, comparison of the 
calculated results was made with data obtained from microwave 


backscattering from deciduous trees. The calculated results 
agreed with the experiment data (in level) . The most inter- 
esting fixture was, that the depolarized component of back- 
scatter was accounted for correctly. This was due to the 
fact that the anisotropic nature of the scattering medium 
was correctly accounted for. 

Finally, we would like to acknowledge the support and 
encouragement of Dr. D. LeVine and Mr. R. Meneghini of NASA. 
We would also like to thank. Mr., Selim Seker of GWU for per- 
forming the numerical calculations. 


INTRODUCTION 


This paper studies electromagnetic backscattering from 
a sparse distribution of discrete lossy dielectric scatterers 
occupying a region V. The scatterers are assumed to have 
random position and orientation. Scattered fields are calcu- 
lated by first finding the mean field and then by using it to 
define an equivalent medium within the volume V. The scatter- 
ers are then viewed as being embedded in the equivalent medium; 
the distorted Born approximation is then used to find the scat- 
tered fields. This technique represents an improvement over 
the standard Born approximation since it takes into account 
the attenuation of the incident and scattered waves in the 
equivalent medium. 

In the past, electromagnetic scattering from a collection 
of discrete scatterers has been modeled by continuous and 
discrete random medium techniques. In the continuous case, 
the random medium is modeled by assuming that its permittivity 
e(x) is a random process whose moments are known. The average 
backscattering cross section is then calculated from a know- 
ledge of the statistics of e(x) . Usually it is only the mean 
and correlation of the permittivity that are required. The 
analysis of this problem can then proceed in a number of ways. 
One method involves calculating the mean field [Keller, 1962, 
Tatarskii and Gertsenshtein , 1963 and Keller and Karal 1966], 
using it to define an equivalent medium and then employing 
single scattering in the equivalent medium. In active remote 
sensing applications from terrain, the technique has been ap- 
plied by a number of authors [Rosenbaun and Bowles, 1974; 
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Hevenor, 1976; Fung and Fung, 1977; Fung and Ulaby, 1978; 

Fung, 1979; and Zuni,ga, et al., 1979]. Another technique 
used to obtain the scattered fields from a continuous random 
medium is the radiative transport approach. Here the trans- 
port equations are obtained in terms of the statistics of e(x) 
[Tsang and Kong, 1978], 

In this paper, we have adopted the alternative approach 
- modeling by discrete random media techniques. Here, the 
individual objects - such as leaves - are characterized by 
their scattering cross sections or dipole moments. Each object 
is then given a random placement and orientation. Techniques 
such as Born approximation; single scattering method and trans- 
port equations have been used to calculate scattered fields 
[Ishimaru, 1978]. In the area of active remote sensing from 
terrain Du and Peake [1969] have used the Born approximation 
to calculate the scattering from a layer of leaves. Lang [1979] 
employed the distorted Born approximation to calculate the 
backscatter from a half space of spherical scatterers. 

In this paper we generalize the previous result of Lang 
[1979] to scatterers of arbitrary shape. The technique employed, 
as mentioned previously, is to find the mean field using a 
technique first developed by Foldy [1945] and later generalized 
by Lax [1951, 1952], Twersky [1962, 1978] and Keller [1964]. 

From this mean wave an equivalent medium is derived which, in 
general, is inhomogeneous, spacially dispersive and anisotropic. 
The scattered field is then obtained by employing single scat- 
tering in the equivalent medium. 

The method is then applied to a half space of scatterers 
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ft 


that are homogeneously distributed and have characteristic 
dimensions small compared to a wavelength. In this case the 
equivalent dielectric tensor is homogenous and nondispersive 
but is still anisotropic. Further simplifications are obtained 
when the orientation of the scatterers is assumed to be dis- 
tributed uniformily in the azimuth direction. In this special 
case the medium becomes uniaxial. By using this equivalent 
medium, simple expressions are obtained for the horizontal, 
vertical and cross polarized backscattering coefficients. 
Finally, the method is used to model a leaf canopy when the 
leaves are modeled by lossy dielectric discs. 
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PROBLEM FORMULATION 




Consider the problem of scattering of time harmonic 

electromagnetic waves from N discrete scatterers located in 

a volume V as is shown in Figure 1. The particles are all 

identical and each has volume V , relative dielectric constant 

P 

e r and free space permeability It is assumed that the 

background medium is free space having permittivity and 


permeability Vq. 


. th 


The position of the i particle is specified by the 

vector extending from an origin 0 to the center of that 

particle c The particle's center is located by the center of 

the smallest circumscribed sphere in which the particle can 

be placed. Although the particles are identical they have a 

rotation with respect to a fixed direction. The rotation for 
th 

the i particle is specified by where 0^ and 

are polar and azimuth angles respectively with O<0^_<7r and 

0_£$i_<2 TT . 

The electric field obeys the vector wave equation 


VxVxE - k 0 e r (x)E = iu>y 0 J 


( 1 ) 


where a time dependence e has been assumed. In (1) k R = 

is the free space wavenumber and J is the current 
density of the source. The relative dielectric constant 
e r (x) can be expressed in terms of individual particles by 
employing translations and rotations of the particle located 
at the origin. Let us assume that a particle located at the 
origin is characterized by the function U (x) where 
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( 1 , xeV 

U(x) =1 P (2) 

( 0 , x£V p 

Using (2) we express £ r (x) as 

N 

e (x) = 1 + A £ U(x-x. ,fi. ) , A = e -1 (3) 

i=l 11 r 

where 

U (x, fi) = u (B(fi)*x) (4) 

Here U(x,fl) is the function U(x) rotated by ft and .g(ft) is a 
rotation dyadic. 

We will find it convenient to express (1) and (3) in a 
more abstract notation. We have 

N 

(L - £ v. ) -E = g (5) 

i=l - 1 - 

where 

L = VxVxI - JCqI (6) 

V i = kQAU(x-X i ,fi i ) I , g = iuy 0 J (7) 

Here I is the unit dyadic and g can be viewed as a normaliz 
ed source term. At times it will be convenient to write 

E(x) = E q (x) + E s (x) (8) 

where Eq(x) is the solution to (5) when no scatterers are 
present, i . e . , 

L-E Q (x) = g (9) 

and E g (x) is the scattered field from the particles. 


5 




. ■■ 


SINGLE SCATTERER - TRANSITION OPERATOR 


Before considering the N particle scattering problem, we 
will consider scattering from one particle located at the 
origin. Putting N=1 in (5) with X^O and we have 

(L - V) -e = , V = AkgU(x,f2) (10) 

where 

e = e Q + e g , t*e 0 = 2 ( n > 

and e is outgoing as |x|->°°. We have used the small e nota- 
s 

tion for the field here to remind us that there is only one 
scatterer present. 

If we use (11) in (10) , we obtain 

= Y*e (12) 


From (12) we see that the term on the left, V*e, can be view- 
ed as the source of the scattered field. We write 

g = V-e (13) 

— eq = — 


where 


g e ^ is an equivalent source term. 

the sources, g , exist inside the 
— eq 


Since Y=0 when 
particle boundaries. 


It is more natural to think of the equivalent sources as 
being caused by the incident field e Q . Because Maxwell's 
equations are linear, we can write 


Seq = 1-So 


(14) 


where the dyadic operator T is known as the transition opera- 
tor in the scattering literature [Lax, 1951]. Now using (13) 
and (14) in (12) and multiplying through by L \ we have 
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( 15 ) 



Thus the knowledge of T completely characterizes the scatter- 
ing properties of the particle. The operator T is related to 
the dyadic scattering amplitude of the particle and for dipole 
scatterers T can be determined from the polarizability of the 
particle. Thus T is directly connected with quantities of 
physical interest. 

The transition operator is a linear bounded operator and, 
as a result, can be expressed in integral form: 

3 e q (x) = l’£ 0 = ) - e Q ( x ' ) (16) 

where the limits for the integral extend over all space. One 
can show that t is 0 wheri x and x' are outside the particle 
[Frisch, 1968], i.e., 

i(x,x') = 0 , x ^ V p or xVV (17) 

The property follows directly from the fact that the equiva- 
lent sources for the scattered field are located within the 
particle boundaries,, 

We will now represent £ in terms of plane waves. The 
representation of t can be directly related to the dyadic 
scattering amplitude. We proceed by representing e Q (x) by 
its Fourier transform, putting this in (16) and taking the 
Fourier transform of (16) . We obtain 

3 (k) =f dk't(k,k') •e 0 (k') (18) 

where 

t(k,k') = — /dx dx , t(x,x , )e" l( --"- , - l) (19) 

(2tt) J J 

In (18) we have used the notation that h is the Fourier 
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transform of h_. More specifically 

h (k) = fax h(x)e‘^- (20) 

Inverting (19) the transition kernel t can be expressed in 
terms of its plane wave representation t : 

t(x,x') = — — — sr /dk dk' t (k,k*)e i( ^*5"^' (21) 

( 2 TT ) J J * 

The dyadic scattering amplitude will now oo. defined. 
Consider a plane wave incident upon a scatterer located at 
the origin. An arbitrary pdane wave can be decomposed into 
two mutually orthogonal linearly polarized plane waves. The 
polarization directions are taken as a p and 3_° where a° and 
3_° are orthogonal unit vectors with a° and 3_° being perpen- 
dicular to the direction of propagation. The two incident 
waves are 

””^ 0 — ' — 

e 0 (x,i;q) = g°e , qe{a,3) (22) 

where i is a unit vector in the direction of incidence. It 
is more convenient to consider both polarizations simultan- 
eously so we introduce the dyadic incident wave [Twersky, 1967] 
e 0 (x,i) = e^ (x,i_;ot) < 23 ) 




w 


ik 0 i*x 


(24) 


iLi*x 

* (I - ii) e U (25) 

The dyadic scattered field from the particle is given by 

i s (x,i) = e g (x,i.;a) a Q + e s (x,i;$)@ 0 (26) 

where e s (x,i;q) is the scattered field due to polarization q. 
The dyadic scattering amplitude, f, is defined in terms of the 
asymptotic expression for § g in the radiation zone. We have 
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( 27 ) 


ik 0 l*! 

e g (x i) <u |(0,i) | — | — , |x|-w» 

where 0 is a unit vector in the x direction, £=x/|x|. 

The relationship between £ and t can be found by employ- 
ing (15) for large |xj (Appendix A) . The result is 

f(£,i) = 27T 2 (i-0 0) •§ (k 0 0,k Q i) • (I-i i) (28) 

From this relation, we see 

0 * f = ° , £*i = 0 (29) 

Thus | is a four component tensor - all combinations of twc 
incident polarizations and two scattered polarizations. We 
also note that £ does not completely determine t but only 
partially specifies it. In particular, a knowledge of £ for 
a free space wave number k Q only determines E(k,k') at jk| = 

|k' j — J -Q ; also only four of the nine components of t in the 
polarization directions are determined. 

Before concluding this section, transition operators for 
particles not located at the origin will be needed. As before, 
the equivalent sources 3 .^ for a particle located at can 
be related to the incident field. It follows that 

■ 2 eq^-^ = =i*-0 = f=i. ‘-0 } d -' ( 3 °) 

By shifting the sources and the incident field to the origin 
can be related to £. One finds 

tf (Xf x‘ ) = £ (x-Xj. (31) 

Note that throughout the discussion the dependence of £ on 
rotations has been suppressed for convenience of notation. 
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COHERENT FIELD 


In this section we will develop an approximate equation 
for the coherent field by employing the Foldy approximation 
[Foldy, 1945] . The equation is in terms of the transition 
operator and thus, when the scattering amplitude is known, 
the equation is completely specified. After the equation has 
been derived it is pointed out that outside V the coherent 
fields obeys Maxwells* equations with free space permittivity 
and permeability. Inside V, the coherent fields obey Maxwells 
equations with free space permeability and a macroscopic per- 
mittivity that is inhomogeneous, anisotropic and spatially 
dispersive . 

Before discussing the coherent field, the statistics 
that govern the particles position and rotation must be specified. It 
will be assumed that the position vectors X^ f i=l...N and ro- 
tation vectors , i=l. .N are random variables that are spe- 
cified by a 5 n dimensional distribution function. In addi- 
tion, it is assumed that interchanging particles leaves the 
distribution function unaffected. From this general distri- 
bution function we can obtain the probability density func- 
» th 

txon for the i particle. It xs 

p X.fl.(-,-) = p xQ*~'-* i=1 N ( 32 > 

where w=(0,4>). In (32) we have explicitly noted the fact 
that the particles are identically distributed by omitting 
the index i on the left hand side of (32) . We will assume 
that the particles location and rotation are independent, thus 

= P x (x)p fl U) (33) 
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with the usual property: 


/v P * ! - 


(x) dx = 1 


f 


p r (w ' dm 

a go “* 


The particle density is defined by 


so that 


p(x) = N p x (x) 


/ p (x) dx = N . 
J v 


= 1 (34) 


(35) 


(36) 


In addition to the one particle density - when treating the 
correlation of the field - the two particle density will be 
required. We have 

A A A A A A 

Px.x.a.B. f 32!'“i'S>) - p x J b £ <£'*:'“■“> = P xfl <£.“)PS$<5-S> (37) 

- l — i — — 

i, j = 1, . . .N 
til til 

In (37) we have assumed that the i and j particles are 
independent. The independence assumption is valid when the 
particles are sparsely distributed; the case we intend to 
treat. 

We will now develop the approximate equation for the 
coherent field. We start by noting that the total field E 
can be thought of as a sum of the incident field Eq plus a 
sum of the fields scattered from each particle, . We 

have 


E = E„ + E E 1 (38) 

- -° i=l 

^ll 

The total field incident on the i particle is called the 
effective field and is denoted by E ^ . Thus T^*E^ repre- 
sents the equivalent sources generated by the incident field 
in the i th particle and 
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( 39 ) 


E*^ = L _ ^ • T . *E ^ 

— s = =x — 

Using (39) in (38) , we have 

N , . . . 

E = E n + E L -T. »E ' 1 
- -u i=1 - -i - 

This is the equation that we wished to obtain. 

Now we average this equation. The result is 

N , , . . 

<E> = E„ + E L •<T.*E Ui > 
u i=l " ~ a 


(40) 


(41) 


To obtain an approximate equation for the mean we follow 
[Foldy, 1945] and assume 

gU) ~ <£> (42) 

This means that the random quantity E^"^ is to first order 
equal to a deterministic quantity, i.e., to first order it 
is a ergodic quantity. Using (42) in (41) and noting that 
< li >~<g^ *<E>>=<1\> *<E> we have the approximate equation 

for the mean field 

N -1 

<E> = E n + E L *<T . > <E> (43) 

u ±=1 = “I 


Denoting explicitly the dependence of upon and av- 

eraging and then using (3 3) , we have 

<% L > = < l(5. i ,fi i )> = f ds dw P XQ (s,w)I(s,w) 



P x (s)J (s) 


(44) 


where 


T(s) 


= f dw 
4 ?f 


p fi (w)T(s,w) 


(45) 


In (45) the bar over T has been used to indicate an average of 
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angular varies only. By putting (44) in (43) , by noting that 
the scattered terras are identical and by introducing p (s) via 
(35) , we obtain 

<E> = E_ + / ds p ( s) L ^*T(s) *<E> (46) 

-'v ----- - 

Multiplying from the left by £ and using (9) , we get 

«C*<E>= g (47) 


where 





p (s) 2?(s) 


(48) 


This is the equation for the coherent field. 

The arguments that have led to the approximate equation 
(47) have been largely heuristic. The essential approximation 
is contained in (42) where the effective field is assumed 
approximately equal to the average field. Although we will 
not discuss the conditions under which (42) is valid, it will 
be shown elsewhere that the approximation is valid when the 
fraction of volume occupied by the particles is small com- 
pared to the total volume, i.e., NV ^/V<<1. We shall refer to 
a distribution of scatterers satisfying this condition as a 
sparse "distribution. 


Before proceeding we will write the equation for the mean 
in more concrete form. Using (30) and (31) in (47) and (48) , 
we obtain 


L*<E(x) > 



p (s) t (>S - s , x ' -s) *< E ( x ' ) > = g 


(49) 


where 


i(x,x') = f 
4tt 


dw 


(w) t ( x , x 1 ;to) 


(50) 
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Here the kernel ^(X/X'jw) is the same as given in (31), how- 
ever we have explicitly shown its dependence on the angular 
coordinate w. 

We can now use (49) to obtain a macroscopic form of 
Maxwell's equations. First averaging the Faradary's law 
equation; we have 

Vx<E(x)> = icojjQ<H(x)> . (51) 

Then by using (51) in (49), we obtain the macroscopic Ampere's 
law equation. 

Vx<H(x)> = J - iw<D> , <D> = £Qe.* < E> (52) 

when e is a macroscopic permittivity operator which describes 
the average behavior of the medium and J=g/(iu)y 0 ) . It is 

e = I + / ds /dx' p (s) t (x-s ,x' -s) . (53) 

= - k 0 2 - J - - ■ 

This expression simplifies to I (free space) when x£V. To 
see this we note that when x^V, we have x-s/Vp si nce seV. 

Now using (17) we have £=0, When xeV (47) does not simply 
in general. It describes an anisotropic, inhomogeneous, spa- 
tially dispersion medium. 

Let us examine how (53) reduces to some more familiar 
expressions in some special situations. We will assume that 
V is infinite through the remainder of this section. First 
we will consider the case when the density is constant, i.e., 
p(i 3 )=p. In this case the permittivity is translationally in- 
variant or homogeneous. To see this we substitute (21) into 
( 53 ) and we perform the integrations over s^ and k' . We ob- 
tain 
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(54) 


£ “ I + 


f&x' y*dk e 1 -' g(k,k) 

k 0 


when §(k,k)= f d w P n ( u ) t (k ,k ' ;u> ) . 

~ - " 


(55) 


Since the integrand is a function of x-x" the permittivity is 
translationally invariant however it is still anisotropic and 
spatially dispersive. 

Another special case of interest is when t is scalar, 
i.e. t(x,x') = t(x,x')2. This occurs when the scatterers are 
spherical. Then the permittivity is isotropic but inhomo- 
geneous and spatially dispersive. 

The last special case to be treated is when the wave- 
length is large compared to the size of a scatterer. Here 
the particle can be treated as an electric dipole. Its equi- 
valent source distribution is given by 

-Seq = iwy 0^eq = iwy 0 (-iw£<$ (x) ) (56) 

where p is the electric dipole ncnient of the scatterer and 6 (x) is tne 
Dirac delta function. The dipole moment is related to the incident field e^ by 
the polarizability tensor g [Jones, 196^ ] 

P = £ 0 =‘-0 (57) 


Using (57) in (56) and comparing it with (14) , we find 


T = k^6 (x) ( 58 ) 

or 

t(x,x') = kQ^6 (x) 6 (x ' ) (59) 

Now putting (59) in (19) , we obtain 

£(k,k') = k 2 Q g/(2ir) 3 (60) 
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Thus we see in the dipole limit t is independent of k and k_* . 

Since we have an expression for % in the low frequency 
or dipole limit, the special form of g can be easily obtained. 
Using (59) in (54), we find 

e = I + p (x) a (61 ) 

Thus in the low frequency case, the permittivity is no longer 
spatially dispersive however it is still anisotropic and in- 
homogeneous . 


CORRELATION 


In this section we will calculate the correlation of the 
electric field. Rather than following procedures used to find 
the coherent wave, the distorted Born approximation will be 
employed. This is a single scattering approximation where 
the scatterers are assumed to be embedded in the equivalent 
medium which has been found in the previous section. The 
method is useful when the fractional volume is small (NVp/V<<l) 
and the albedo of a single particle is small. The later con- 
dition implies that the energy absorbed by a particle must be 
much larger than the energy scattered by it. 

We start by considering a volume V of equivalent medium 
surrounded by free space. There are N particles embedded in 
V as shown in Figure 1. The scattered field due to the i^ 
particle can be calculated by modifing (39) . We assume that 
the incident field on the particle is the mean field <E> and 
that the free space operator £ is replaced by the equivalent 
medium operator 3? as given in (48) . We have 

N , . , N , 

E = E E u ' = E • (62) 

i=l ~ s i=l " =1 - 

Before proceeding, we point out that our main interest 
in finding the correlation of the field is to use it to cal- 
culate the backscatter ing cross section. Since this cross 
section is related to the correlation of the field fluctua- 
tions, we now define 

E f = E s - <E s > , <E f > = 0 (6 3) 

Now computing the correlation of the fluctuating field, we 
obtain 
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(64) 


< E f (x)E*(x) > = <I s (x)E*(x)> - <E g (x) ><E* (x) > 

where z* is the conjugate of z. Putting (62) in (64) and no- 
ting that a portion of <E.E*> cancels with <E ><E*> if we use 

-S -s — s — s 

the fact that N>>1. We find 

< E f ( x) E * (x) > = <% P n (o))<E f (x)E* (x)> w (65) 


where 


E f (x)E*(x) = J ds p(s) £ (x,s)£* (x,s) ( 66 > 


with 


3-1. 


(67) 


£(x,s) = f *T(s) -<E> 

Here we have separated the average into rotation and coordi- 
nate space averages, thus introducing the conditional expecta- 
tion, <E f E|>^, with respect to given w. 

To write (67) more explicitly, we introduce the dyadic 
Green's function g(x,x') for the operator It satisfies 

f-g(x,x') = |6(x-x') (68) 

+ g - outgoing as |x|-*<» 
where £ is given in (48) . Now (67) becomes 

£(x,sj = fdx'g{x,x') • ydx"£(x , -s,^'-s) -<E(x")> (69) 

The expression simplifies greatly in the low frequency limit. 
Assuming that t is given by (59) and using this in (69) gives 


£(x,s) = (x,£) *^*<E (s) > 


(70) 
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BACKSCATTERING COEFFICIENTS FOR A HALF SPACE OF DIPOLES 

To illustrate the application of the methods developed 
in the previous sections, we will calculate the backscattering 
coefficients from a half space of scatterers that are small 
compared to wavelength. We will also assume that the density 
of scatterers p is constant. The physical configuration is 
shown in Figure 2. There, we have shown the direction of the 
incident wave and the polarization vectors h° and v° repre- 
senting horizontal and vertical polarizations respectively. 

To compute the scattered field using the distorted Born 
approximation, we must first calculate the mean field in the 
half space containing the particles. In the low frequency 
approximation the mean wave is computed by replacing the par- 
ticles with an equivalent medium having relative permittivity 
tensor g=|+p g and free space permeability . The usual con- 
tinuity conditions associated with macroscopic Maxwell's 
equations are assumed to hold at the interface z=0. 

Before proceeding we would like to emphasize that the 

scatterers are sparsely distributed or that the fractional 

volume they occupy is small (NV /V<<1) - a condition necessa- 

]? 

sary for the validity of the mean equation. This restriction 
is reflected in the equivalent permittivity tensor. Small 
fractional volume requires that |poUj|<<l where are the 

components of g. We can exhibit the dependence of this con- 
dition on the fractional volume NV p/ v= P V p explicitly by in- 
troducing a normalized polarizability tensor g as follows: 

a = g/V p (71) 
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where we can show that the components of a remain bounded as 
V -*0. Now the permittivity can be written as 

fr 


| ^ I + 6a (72) 

and thus we have a small parameter for ordering purposes . 

Although we are able to carry out the calcualtion of the 
mean wave for an arbitrary average polarizability tensor a, 
it is convenient to partially specify the angular probability 
density, p^(w) in order to make a diagonal. First we choose 
a spherical coordinate system of mutually orthogonal unit 
vectors r Q ' io and . The position of these vectors is com- 
pletely determined by the spherical angles 6 and $ as shown 
in Figure 3. Now we align the principal axes of the scatter- 
er along these unit vectors. Then we write 


§ - a r r°r° + a 0 _e o 6_° + a.^°^° (7 3) 

By using the usual transformation between spherical and car- 
tesian coordinates, (7 3) becomes 

3 3 

a - L E a v (9,(J>)x?x? (74) 

i=l j=l x i x j 

where x 1 =x, x 2 =y and x 3 =z and x°,v° and z° are cartesian unit 
vectors. The relationship between the a x x and a r ,a 0 ,a^ are given in Appendix 
Now assuming that the random variables 0^ and are in- 
dependent, we write 

P„<“> = P 0 ( £l )p t ( , H ' (75) 


Averaging (74), we have 


a 


3 3 

E E 
i=l j=l 


a x?x° 


(76) 


where 
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a x x = J M J d'J) p q( 0) p .(<})) a (Q,$) (77) 

x i i 00 9 x i x j 

We will now assume that the particles are distributed 

uniformily in the <j> variable, i.e., p* (<j>) =1/2 tt . We find 


a = a 


xx- — + a yyi“i. 0 + a zz Z°*° 


where 


a xx = a yy ' l [a r si " 20 + a 6 + V (79 > 

a zz = a r cos2e + a e sin2 9 (80) 

Thus all off diagonal terms averaged to zero and two of the 
on diagonal terms are equal. Using this in (72) , we see that 
the equivalent dielectric is uniaxial. 

The mean wave in the equivalent medium will now be found 
for the case of particles uniformly distributed in the azimuth 
coordinate cj> . The incident wave is given 

E Q (x,q) = 2 °e 1 -*- qe(h,v} (81) 


where 


with 


k = k, + k z‘ 


\ = k o sin0 o x" , = k 0 cos6 0 


and the polarization vectors are 

h° = , v° = -cosBqX 0 + sine^z 0 (8 4) 

The average electric field in the equivalent half space satis- 


fies 


[Vx(VxI) - (X+6a) ) *<E (x) > = 0 , z< 0 


Let us assume a plane wave solution of the form 
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<E (x) > = A e 1 -'- (86) 

where k = k .. + K z° . In order to match fields at the inter- 
— — t z— 

face, the transverse phase velocity of the incident and trans- 
mitted waves must be the same, thus ic_ t =k t =kQsin0QX° . Putting 
(86) in (85), we have 

k_x (kxA) + kgA + 6k^a’A =0 (87) 

Representing A in cartesian components (87) can be written as 


K z‘ k 0 B xx 


k 0 sin8 0 K z 


k O sln0 O K z 


2 , 2 „ 

K z*Vyy 


- k 0 B zz 


0 (88 ) 


where 


8 = 1 + 6a , 

XX XX 


0 = cos 8 n + 6a , 

yy o yy 


8„_ = COS 6 n + 

zz 0 zz 


Since (88) is a homogeneous system, the determinant of coef- 
ficients must be zero for a solution to exist. This condition 

determines the allowable values of k_. We find 

z 

* z = ±4 M - ± k </v <90) 


- ±<4 V> = ±J- 

.si 


0 

XX zz 


sin V 3 zz- 


where the superscripts h and v have been used to designate the 
propagation constants associated with horizontal and vertical 
polarizations. For an incident wave that is not grazing, i.e. 
SqA/ 2, expressions (90) and (91) can be simplified using the 
small 6 parameter. We have 
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(h) _ 


6 a 


= k~ (cose. 


l 0\ wou 0 + 2 cos 07/ + e(6 " ) 


(92) 


(v) _ 


[= 


s> /a +cos^ 8 n (a -a )) 

k Icosa + °{ 22 0 xx zz 'l 

k rt |c°s0 o + ooiej ) 


+ 0(6*) (93) 


We cannot calculate these propagation constants to higher ac- 
curacies than 0(5) since the original mean equation has only 
been found to this accuracy. We note that since we are con- 
sidering lossy particles, the are complex and thus the 

z 

mean wave will decay away from the interface. 

Next we calculate the amplitude coefficients for the 
mean wave of both polarizations. We have 

ik z z -ik 2 z — 

(e 6 +r e 0 ) e y° , .z <0 

h 

(94) 

• (h) 

1X_ V ' -X 

A e y° , z>0 

y 


<E(x,h)> = 


and 


ik z z -iJc z ^t *» 

(e ° +r e ft ) e v° , z <0 

V (95) 

• ( v) 

(A x x°+A z z°) e 1 - , z>0 

where we have introduced a reflected wave in the free space 
medium at the specular angle. Now by using the fact that the 
tangential <E> and <H> must be continuous at the interface, 
the unknown reflection and transmission coefficients can be 
calculated. Since the major effect of the equivalent medium 
is to produce exponential decay, we next expand the coeffi- 
cient for small 6 and we keep only zero order terms. We find 
that r = 0 ( 5 ) and thus it can be neglected. The transmitted 

Si 


<E (x,v) > = 
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jjasBBKBg SiS a 




mean fields are: 

i(K z ^z+k fc -x) 

< E(x,q)> = c[°e 0 + 0(6) , z>0 (96) 

qe {h, v} 

Proceeding with our developemnt, we now relate the transverse 
Fourier transform of the correlation to the backscattering 
coefficients. This is done within the context of the distor- 
ted Born approximation developed in the previous section. We 
start by taking the transverse Fourier transform with respect 

A 

to x and x of (65), (66) and (70). We have 


cEf (k t ,z)E*(k t , 2 ) > = 

*'4 ir — — 

(97) 

<If (]St ,z) -f ( ^t ,z) > U = 

/• -v ^ A A 

p / ds £(k. , z,s) •£* (k., z,s) 
* / V z 

(98) 

where 




£(k. t fZ,s) = kQgtk^ZitS) *a*<E(s) > (99) 

Since we will only require Q when x and s are in the equivalent 
medium and since the reflection at the interface is small we 
can replace g by the dyadic Green's function for an infinite 
equivalent medium, i.e., 

S(x, s) = S (K>) (x-s) + 0(6) (100) 

We have written the infinite space Green's function in terms 

of x-s since it is translationally invariant. We then have 

— ’ / *** j \ — ik, * s 

i (k t ,z,s) = 9 {m) (k t ,z-s)e ~ (101) 

Now by putting (9g) , (Ss) , (100)and (101)in (98), by integra- 

A 

ting over s fc and by setting z=z=0, we have 

< E f (k t »0)E*(k t ,0)> w = i(k t ,q|u))6(k t -k t ) (102) 

where 
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g(k t ,q|u)) « ( 2 TT) 2 £!^V p 


/: 


ds (g (oo) (k fc f -s) '§*a°) (S ( °^ * (k fc / -s) 


-2ImK (q) s 

*§*'S °) -e 2 (103) 

Here §(k. t »qjw) is the transverse dyadic spectral density at 
the interface assuming w is fixed. The normalized polariza- 
bility has been introduced by using ( 71 ). 

By using the results of Appendix C and by noting that 
ITk^Tq^jT we obtain the backscattering coefficients 


cr 


k o c ° s2(5 o 


pq 


— 3 — £°*§(“li t /q|a)) *P° f p # qe {h,v3 
4 tt 0 


(104) 


To evaluate the integral of £ in (103) , we will need the 
transformed Green's function. To obtain it, we first write 


the governing equation for G 


(«) 


It is 


[Vx(VxI) - k 2 (I+<Sa)] * G (00) (x) = | 6 (x) , 


Using 


, (°°) 


,(°°) 


(x) - outgoing as lx )-*- 00 


(2tt) 


f 


(X) = — ^ JdK 


in (10 5), we find 


(105) 


(106) 


[(kxjocI) + k 2 (I+6|)]*g(^) = -I 


Then 


- h 


h 


die g(K)e^^ 


(co) 


(107 ) 


(108 ) 


To simplify the remaining computation for Q we note from 

(00) 

(104) and (103) that G will only be required for £=-k. = 

t 0 

-k o sin0 o x° . 

Inverting (107) and performing the integral in (107) by 
the method of residues, we obtain 
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S M (-K. .*) 
r o 


= 1 


sin0 n K^ v) o (z) 

S x°x° + (x°z°+zV) 

ZZ K ~ — — — — 


‘* 4 ” - k 0 6 x X ) 


) 2iK 




(v) 


iK lh) |Z| 

+ fi>> 


( 109 ) 


where 


cr(z) = 


( 110 ) 


j 1 , z>0 

(-1 , z< 0 

and the 0's are defined in (89). Approximating the coeffi- 
cients to zeroth order in 6, we have the simplified expression 


, i 4 V) I z I 

1 ( -k , ,z) = V°v°^ -- 

-t Q - - 2l k o cos0 o 




+ h°h° 


— 2ik o cos0 Q 


(111) 


If we use this in (103) , perform the integration and use 
the result in (104) , we have our final form for the backscat- 
tering coefficients. It is 


o° 

pq 


5k» V la | 

2. pj-pa. 1 


8i7(Imk^ p ^ + ImK^) 
z z 


( 112 ) 


where 


I a D a 1 2 = f dw|£°*a-2°| 2 . (113) 

pq 

and < v , se {h ,v} are given in (92) and (9 3 ). The dependence 
z 

2 

of I a pg I on angle of incidence is worked out explicitly in 
Appendix B for scatterers that are distributed uniformly in 
the <j) coordinate. 

The final result given in (112) can be expressed in terms 
of scattering cross sections of individual particles. Using 
(28), (60) and (71), we have 
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C3 


f „ = p° • f *q c 

pq * = M 


2tt 2 S = k^V a /4 t t 
pq Op pq' 


(114) 


Now by recalling that the backscattering cross section from 
a particle a' -4ir|f I , (112) becomes 

pq 1 pq 1 


pa (b) 

a o = 7— 


2Im« +2lmK 
z z 


p,qe {h,v} 


(115) 


Following Attema and Ulaby [1978] we can give a one dimen- 
sional interpretation of (115) . If we rewrite (115) as 

f 0 -2lmK^ ul 

ct° = / dz pa W) e a 

pq M pq 


z | -2lmK i ; i ' | z 
e Z 


(116) 


we can view the scattering as being decomposed into scattering 
from slabs of width dz. An intensity of exp(«-2ImK^ |z|) is 
incident on the slab located at z. The incident intensity is 

/ i_ \ 

backscattered with reflectivity factor • The backscat- 

tered wave then decays as exp(-2Imt<^ |z|) until it reaches 
the interface. 


27 


DISCUSSION AND NUMERICAL EVALUATION 


In this section, we will first discuss several general 

properties of c^ that are independent of the particular scat- 

terer chosen. Following this discussion, we use our method 

to model a forest canopy by a collection of lossy dielectric 

discs. The theoretical curves computed from this model are 

then compared with some experimental data. 

Because of the simple dependence of a° on the medium 

pq 

properties and incidence angle certain general observations 

can be made that are independent of the particular nature of 

the scatterer. First, we note that ff°g as given by (115) is 

independent of the density of scatterers p. This follows di- 

( s ) 

rectly from (92) and (93) where we see that the Pk , se {h,v> 

z 

are directly proportional to p. Thus the linear p dependence 

in the numerator of (115) is cancelled out by t-he denominator. 

Second, we note that " cos8 ft . This is the same result as 

hh U 

predicted by the scalar theory. Finally, we note that 0° = =a^. v 
at normal incidence (0°=O) . This is an expected result. Since 
the scatterers are uniformly distributed in <p , the two polari- 
zations see the same medium at normal incidence. 

We now proceed to model a forest canopy by a collection 
of leaves. The leaves are in turn assumed to be lossy dielec- 
tric discs as mentioned previously. The discs have radius a 
and thickness h. Typical dimensions are radii of one to sev- 
eral centimeters and thicknesses of tenths of a millimeter. 

The electrical properties of discs can be characterized by 
their normalized polarizability tensor a when the wavelength 
is large compared to the disc. From Jones [1964] and 
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Van de Hulst [1957] the polarizability of a disc along its 
principal axes is given by 



A 

1+A ' 


a e 


= A , 



(117) 


when r ,6,<p are defined in Figure 3. 

Because of the large volume of water present in vegeta- 
tion, we can usually assume | | > > 1 in the microwave region. 

Using this assumption in (117), we find that [ a^ |= j a^ | >> | | . 
This inequality can now be used to simplify the scattering 

cross section of (115). We find that a 0 «:|e J 2 v /Im£ . 

PtI P ^ 

Thus it follows that the magnitude of the backscattering cross 
sections are directly related to the volume and complex die- 
lectric constant of the discs in a simple manner. Therefore, 
as leaves grow and as their moisture content changes these 

effects should be observable by measuring a 0 at different 

Pm 

periods of the growing season. 

Before computing the backscattering cross sections as a 
function of incidence angle, we will require the relative di- 
electric constant of the leaves and the angular distribution 
of leaves. First the relative dielectric constant is considered 
Our calculation follows that of Fung and Ulaby [1978] who in turn have 
based his results upon de Loor [1968] and Carlson [1967]. 

They model the leaves as a mixture of water and solid materi- 
als. For illustrative purposes we have chosen 50% water and 
50% solid for our calculations. By using (3) and (4) of 
Fung and Ulaby [1978] at a frequency of 1.1 GHz we find that 5^= 
30.8+il.8. Our choice of frequency has been motivated by experimen- 
tal results that appear in the literature. We have chosen to 
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compare our results with those of Bush, et al [1976] who has 

measured a®, p,qe {h,v} from forests for frequencies 1-18 GHz. 
PS1 

Because of the dipole approximation made in our model, only 
the lowest frequency (1.1GHz) Ulaby measured was used for 
comparison purposes. 

The angular distribution of leaves will now be considered. 
Field measurements of leave orientations have been made by 
Smith [1973] and others. It has been found that the leaves 
are distributed uniformly in the $ coordinate (Figure 2) . The 
distribution of leaves with respect to 0 is more vegetation 
type dependent. Several are given by Smith [1977]. Since 
no measurements of this type exist for the Ulaby data, we 
have assumed that $ is uniformly distributed. For 0 we have 
considered the following two density functions: 


P e (6) = 



0< 6<_A0^ 

A0 <_0<_7T 


( 118 ) 


or 


( 2A0T ' !L“A0j.l0il + A0 i 

P Q (0) = ( ( 11 9) 

0 , elsewhere 

In (11 8) when A0^, is small, the leaves are approximately 
parallel to the interface (z=0) ; when AB^ =tt/2, they are un- 
iformly distributed in 0„ In (119) when A0^ is small-, the 
leaves are perpendicular to the interface; when A6^ is tt/2 
they are uniformly distributed in 0. 

The numerical calcualtions are presented in Figures 4-9. 
In these figures the backscattering coefficient is plotted as 
a function of the angle of incidence 0g. In Figures 4-7 we 
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have used the angular distribution given in (118) while in 
Figures 7-9 we have used (119). Figure 7 corresponds to a 
uniform distribution over all 9 and thus, for this case, 

(118) and (119) give the same results. One should note that 
since |e r |>>l, a change in f, or just tuifts the level 

of the curves but it does not change their shape. Their shape 
only depends on Pg(0). 

The following trends are observed in Figures 4-7: First, 

is always greater than u^. Their difference increases 

as A 0 a becomes smaller. Second, the cross polarized back- 

scatter becomes smaller as A0^, becomes smaller. In Figure 

7-9, we observe that: First, becomes greater than as 

AS, is decreased; second, the curve for a° tends to flatten 
1 hv 

out as A 6^ is increased; and third, the difference a^-cr^ at 
8^=0 becomes smaller as AS increases. 

A comparison of our theory with the experimental results 
of Bush, et.al [1976] is made in Figure 10. There, we have 
plotted our Figure 8 along with his data for Kansas deciduous 
trees measured in the springtime at a frequency of 1.1 GHz. 
Figure 8 was chosen since it most clearly appears to follow 
the trends of the data i.e., flat cross polarization and a^. v > 

a hh* 

Although our theory follows the trends of the data, it 
is clear from the results that additions to our model should 
be made. The fact that CT a ^ is most likely due to the 
vertically oriented tree branches other than due to the leaves 
that tend to be parallel to the interface. In addition, an 
examination of the numerical results shows that the skin depth 
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for the mean wave is large. Thus at a frequency of 1 GHz the 
underlying ground should be taken into account. 
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APPENDIX A 


Relationship Between the Transition Operator and the Scatter- 
ing Amplitude 

■V 

To find the relationship between f and t, we start with 
(15). Using (23) and (26), we can write (15) in terms of 
dyadic incident and scattered wave 

I s = (1A) 

Next we use (16) in (1A) along with the free space dyadic 
Green's function £ for L " . We have 

i s (x,i) = /dx'r(x,x')/dx" t (x' ,x") 'i 0 (x 1 " ,i) (2A) 

Here the free space dyadic Green's function is given by 

ik 0 ix-x' i 

£(M'> = <3A) 

0 

To obtain § s in the radiation zone, the far field expression 
for X will be required. It is [Twersky, 1967] 

e ik ° iSl 

r(x,x') % (|-00)e U -rn , |x|+~ (4A) 


Now putting (4A) and (25) in (2A) 


i r XK Q 

e s (x, i.) ^ 4 ^: J dx' dx"X (x ' ,x") e 


ik A (i*x" - 0 *x') 


• (I-ii) — 


ik 0 I- 


x 


(5A) 


Finally employing (19) in (5A) , and comparing with (27) we 
obtain the required result 

|(0,i) = 2tt 2 (1-0 0) *t(k 0 0,k 0 i) • (I-i i) ( 6 A) 
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APPENDIX B - Polarizability Statistics 


In this appendix calculate the mean square polarization 
statistics used in (112) . The calcualtion will be performed 
with the assumption that the scatterers are uniformly distri- 
buted in 

First we derive the components of the polarizability 
tensor in cartesian coordinates in terms of the principal axis 
components. The unit vectors r 0 , _0° and are related to x°, 
Y° and z° as follows: 

r° = sin6cos<J>x° + sin0sin<})^ 0 + cosGjj 0 

0_° = cos0cos(J)X 0 + cos0sincf)y o - sin0^° (IB) 

£ 0 = -sin<j)x° + cos<J>y° 

Using this in ( 73 ), the cartesian components of (74) are: 

2 2 2 2 
a = (a sin 9 + a.cos 0)cos d> + a, sin d> 
xx r 0 T <p T 

2 2 2 2 

= (a r sm 0 + a @ cos 0)sin $ +a^cos <j> 

2 2 

a = a cqs 0 + a sin 0 

z z i y 

(2B) 

2 2 

a„„ = [(a sin 0 + a.cos 0) - a , ]cosd»sind) 
xy it y cp 

a = (a -a Q ) sin6cos8cos4> 
xz r o 

a = (a -a 0 ) sin0cos0sin<j> 
y z t u 

The other components are gotton from the fact that § ‘is a 
symmetric dyadic. 

Now we will obtain the components of a in the polariza- 
tion directions in terms of the cartesian components. We 
have using (8 4) 


a hh = — ° ° = a 


yy 


= h°>a*v° = -cos0 Q a yx + sin0 n a 


0 yz 


(3B) 


dL — JUI- ■ d - 1 

nv — * = - 

a vh - v° *a*h° = a hv 

a vv = = “^Vxx " 2cose O sln9 O a xz + sin2 ®0 a zz 

Following this, the mean square polarizabilities can be 
found. They are 


hh 


yy 


I a hv I 2 = cos2e el a yxi 2 + sin2e o! a yzl 2 


‘vh' 


= a 


hv 


2 4 _ | 

a vv = cos 0 n I a 


0 1 xx 


| 2 + sin 4 0 o |a zz | 2 +cos 2 0 o sin 2 0 n (4|a xz | 2 


+ 2R a a* ) 
e xx zz 


where the uniformity in <p has lead to 


a* = a 


n* = 


XX xz 


yy yz 


a* = a 


yx yz 


zz xz 


0 . 


The quantities |a xx | 2 , |a yy | 2 , |a zz | 2 , |a yz | 2 , | a xz | 2 

and a a* can be obtained easily from (2B) by using the 
xx z z 

known statistics of 0 and $. 
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APPENDIX C - Relationship Betvreen the Backscattering 

Coefficients and the Transverse Spectral Density 

We start by considering the fluctuating portion of the 
scattered field, E f , as defined in (63) in the region z_>0 . 
This field can be viewed as arising from sources on the in- 
terface. So that far field quantities can be found, we ini- 
tially consider that portion of E f that arises from sources 
contained within a finite region A on the interface. The ra- 
diated field from the region A will be denoted by E f (x,q) . 

The field E^ (x,q) can be related to the interface fluc- 
A 

tuations by employing a plane wave expansion in the region 
z< 0 . We have 


5 f A-' q> = d — t 


-ik„z+ik t (x t -x^) 


, z<_0 
(1C) 


k z = /k o-IHt| 2 ' lmk z >0 

where E f (0,x t ,q) i s the fluctuating field on the interface 
due to an incident wave of polarization q. The k^ integral 
in (1C) can be asymptotically evaluated for large |x| [Collin 
and Zucker, 1969] . We find 


E (x,q) 'v 


+ik 0 cose 0 + i* Q |x| 
2 tt I x I e 


ik . • x 


f — t — t 

/ dx^ E f (0,x^,q)e 0 


where x has been specialized to the backscatter direction 6 ( 
and k is given by (83) . 

The backscattering coefficients are now defined. They 


0 , . M*!l (P,q) 

cr° = lim lim — — - — 


A-*-oo x ->oo 


AI i (q) 


q / p e { h , v } (3C) 
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where I^(q) is the incident intensity per unit area with po- 
larization q, I (p,q) is the average intensity at the obser- 
vation point with polarization p due to an incident wave with 
polarization q. In view of (81), I i (q)=l. We also have 

I (P#q) = < I E (x,q) * p° | 2 > (4C) 

A 


We can complete the development by first representing 
the field on the interface by its Fourier transform, E f (k t< q) 


E f (0,x t ,q) = — f E f (k|,q)e 


i — t'— t 


(2tt) 


dk' 
— t 


(5C) 


Now by using (2C), (4C) and (5C) in (3C), we have 


pq 


lim lim 


? 2 
k^cos^Q 


A-* 03 x ■+■<» 2 


cos ti n r * 

~5 J dkidk? / dx ' dx "< | £° * E A k ' , q ) 

IT A * / AxA t C f r 


E*(k",q)'£°j >e 




(6C) 


Next we introduce the transverse dyadic spectral density 
S(k t ,q) which is given by 


<E f (k£,q)E*(k£,q)> = § (k^,q) 6 (k^-k") (70 

Note that this definition requires that the process E f (k t ,q) 
be homogeneous or stationary in k fc . 

By using (7C) in (60 and carrying out the integrations, 
we obtain 


where 


pq 


.2 2 C 

V os e o 

— 7~1 

4 TT 


S ( -k 
p -tf 


#q) 


s p(-t' q ) = ■§ { ^t ,q) *E° 


(8C) 


(9C) 
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